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ABSTRACT 



We report on very high energy (E > 100 GeV) gamma-ray observations of 
V407 Cygni, a symbiotic binary that underwent a nova outburst producing 0.1- 
10 GeV gamma rays during 2010 March 10-26. Observations were made with the 
Very Energetic Radiation Imaging Telescope Array System during 2010 March 
19-26 at relatively large zenith angles, due to the position of V407 Cyg. An im- 
proved reconstruction technique for large zenith angle observations is presented 
and used to analyze the data. We do not detect V407 Cygni and place a dif- 
ferential upper limit on the flux at 1.6 TeV of 2.3 x 10 -12 erg cm -2 s _1 (at the 
95% confidence level). When considered jointly with data from Fermi-LAT , this 
result places limits on the acceleration of very high energy particles in the nova. 

Subject headings: (Stars:) white dwarfs — (Stars:) novae, cataclysmic variables — 
Gamma rays: general 



1. Introduction 



In March of 2010, the Fermi-LAT Collaboration announced a new GeV transient in 
the galactic plane, FGL J2102+4542, that was identified as a nova outburst in the sym- 
biotic binary V407 Cygni (hereafter, V407 Cyg). At least seven GeV transients located 
near the Galactic plane have been discovered by EGRET, Fermi-LAT and AGILE. Only 
two have been identified at other wavelengths: V407 Cyg, which is the first no va to be de- 



tected at GeV energies, and J0109+6134, which was likely a background blazar ( lAbdo et al. 
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2010h . The physical nature of the other five sources is unknown 



20101 ; lHavs et all 120101 ; iSabatini et all 120101 ; iTavani et all 11997 



20101 ). and some of these GeV transients may represent a new class of 



gamma-ray emitting objects. 

The Fermi-LAT Collaboration reported variable gamma-ray emission in t he 0.1-10 GeV 



band from FGL J2102+4542 during 2010 March 10-26 (MJD 55265-55281) (ICheung et al. 



20101 ) . Its flux in gamma rays, binned on a day-to-da y basis, peaked 20 10 March 13-14 with 
a flux of 9 x 10~ 7 photons cm" 2 s" 1 above 100 MeV (Ubdo et alJboid ). The GeV gamma- 
ray activity lasted approximately two weeks. The initial report of GeV emission triggered 
Very Energetic Radiation Imaging Telescope Array System (VERITAS) observations of the 
object at very high energy (VHE; E > lOOGeV) as part of an ongoing campaign to observe 
transients detected by Fermi-LAT . 

Using multi-wavelength data, it was determined that the new transient was most likely 
associated with V407 Cyg, a binary system c onsisting of a Mira-type pulsating red giant and 
a white dwarf companion (jAbdo et al.ll2010l). A nova outburst from V407 Cyg was detected 
in the optical waveband on 2010 March 10 (INishiyama fc Kabashimall2010l ) with a magnitude 
of ~ 6.9, while pre-outburst magnitudes from the previous two years of monitoring ranged 
between magnitude 9 and 12 (jAbdo et al.ll2010f ). V407 Cyg has been optically monitored for 
decades and has experienced previous ou tbursts, but the system had never been observed 



to be as bright as during the nova (e.g., iMunari. Margoni. fc Stagnilll990t iKolotilov et al. 



19981 ). The onset of the optical outburst corresponds to the first significant detection of the 
source by the Fermi-LAT on 2010 March 10. 

Novae in red giant/white dwarf systems have been known to pro duce expanding shocks 



that can result i n X-ray emission (e.g., the recurrent nova RS Oph, ISokoloski et al.l 12006 



Bode et al. 



2006h. and indeed, X-ray emission from V407 Cyg was detected after the nova 



( I Abdo et al .1120101 : iNelson et al.ll2012[ ). Based on the observed X-ray emission from the 2006 
nova outburst of RS Oph, before the launch of Fermi-LAT , iTatischeff fc Hernanzl (120071 ) 
suggested that particles could be accelerated in novae up to TeV energies, but gamma-ray 
emission from a nova had never previously been detected. Here, we discuss the VERITAS 
observations of V407 Cyg and their implications for gamma-ray emission from the nova. We 
also describe an improved event reconstruction technique for stereo observations by imaging 
atmospheric Cherenkov telescopes (IACTs) made at large zenith angles (LZA). 
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2. Observation and Analysis 

2.1. VERITAS Observations 

VERITAS is a ground-based VHE gamma-ray observatory located at the Fred Lawrence 
Whipple Observatory in southern Arizona. It consists of four IACTs sensitive from approxi- 
mately 100 GeV to above 30 TeV. Each VERITAS telescope has a 12 m tessellated reflector 
with a total area of 110 m 2 . Each camera's focal plane contains 499 closed-packed circular 
photomultiplier tubes, giving a total field-of-view of 3.5°. Gamma-rays incident onto the 
upper atmosphere induce a particle cascade, called an air shower, in which some charged 
particles have sufficient speed to emit Cherenkov light. The direction and energy of the 
original gamma ray can be reconstructed from images of the Cherenkov light recorded by 
the telescopes. When observing at small zenith angles (< 40°), the array has a n energy reso- 



lution of 15% at 1 TeV and an angular resolution of better than 0.1° at 1 TeV ( lAcciari et al. 



20081 ). For observations at LZA, the energy and angular resolution are degraded and the 



energy threshold is increased. 

VERITAS observed V407 Cyg for several nights after the announcement of the Fermi- 
LAT detection, during days 9-16 of the outburst (2010 March 19-26). The zenith angle 
of these observations ranged between 50° and 66°. The VERITAS telescopes are regularly 
operated in a mode called wobble mode, during which the location of the object to be ob- 
served is offset from the center of the fi eld of view (FoV) b y 0.5°, allowing for simultaneous 



source and background measurements (IFomin et al.lll994l ). The offset direction cycles be 



tween north, south, east and west for sequential observing segments to reduce systematic 
effects. After filtering the data for contamination due to poor weather or instrumental prob- 
lems, 304 minutes of live time remained from the original 335 minutes of observations, see 
Table 1. 

To test the improved reconstruction technique discussed in Section 12. 2\ VERITAS ob- 
servations of the Crab Nebula were also analyzed. We selected 203 minutes of good time 
intervals from 17 data segments taken on the Crab Nebula during 2010 March 12-16 (MJD 
55267-55271) with similar zenith angles ranging from 5 5° to 65°. All dat a were analyzed 



using the standard analysis package for VERITAS data ( jCogan et al.l 120071 ) . 



2.2. Event Reconstruction 



The raw data were calibrated and cleaned, and quality selection criteria based on the 
number of photomultiplier tubes contained in the images and the position of the image in the 
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camera were applied. The s hape and ori entation of the gamma-ray images were parametrized 
by their principal moments (lHilraslll985f ). In order to produce gamma-ray images of the sky, it 
is necessary to reconstruct the putative source location for each shower in the camera plane 
(hereafter "arrival direction"). When imaging showers with multiple IACTs, the arrival 
direction of a shower is usually found using simple geometric arguments. The major axes 
of the images produced by a shower in each IACT camera intersect near the location of the 
arrival direction. The shower arrival direction is calculated by minimizing the perpendicular 
distance to each image's semi-major axis, weighted by the size of each image. This method, 
here called the standard method, is effective at small zenith angles. However, at LZA, the 
major axes of the air shower images from an individual gamma-ray event are generally close 
to parallel. Thus, the uncertainty of the intersecting point increases, resulting in a loss of 
angular resolution. Due to this effect, a reconstruction technique that does not depend on 
the intersection of the axes is desirable for LZA observations. 

The disp lacement method is a direction reconstruction algorithm that is useful for LZA 
observations (jSentiirk et al.ll201ll ). In very general terms, it consists of calculating the arrival 
direction using the shape and brightness of a given air shower image. More specifically, it 
relies on the determination of the disp parameter, defined as the angular distance from the 
image ce ntroid to the arrival direction. This method was used by several e xperiments in 
the past (jLessardll200ll ; iKranich et al.ll2003l ; iDomingo-Santamaria et al.ll2005l ). with varying 
ways of calculating disp. 

The basis of t he displacement meth od is the relationship of the disp parameter to other 
image parameters (IHofmann et al.lll999l ). The implementation of the algorithm in VERITAS 
is as follow s: we estima te disp as a function of three other image parameters, size, length 
and width (iHillasi Il985l ). using Monte Carlo simulated gamma-ray showers. The method 
results in two different arrival directions, one on ea ch side of a telescope image along the 
semi- major axis, also known as head-tail ambiguity (IHofmann et al.lll999l ). This ambiguity 
is eliminated by choosing the cluster of arrival directions closest to one another, one coming 
from each image. Finally, the arrival direction is estimated independently for each telescope 
image and an average weig hted by size is tak en. This method proves to be more powerful 
than the standard method (jCogan et al.ll2007l ) when reconstructing events with zenith angle 
larger than 50°. Quantitatively, an improvement of ~30% in detection significance for a 
source having 1% of the strength of the Crab Nebula has been observed. 
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2.3. Event Selection 



The cosmic ray background rate for IACTs is typically more than 10 4 times the gamma 
ray rate, so it is important to reduce this background while retaining as many gamma-ray 
events as possible. By exploiting the differences in the development of gamma ray and 
cosmic ray induced showers, the background due to cosmic rays can be reduced significantly, 
while still retaining a high fraction of gamma-ray like events. The background reduction is 
performed by placing standard selection criteria, optimized using Monte Carlo simulations 
and real data from the Crab Nebula on the shower image parameters. The selection criteria 
for the size of the telescope images, the mean scaled width and mean scaled length parameters 
( lDaumlll997t iKrawczynski et al.ll2006l ). the height of maximum Cherenkov emission and the 
angular distance from the anticipated source location to the reconstructed arrival direction 
of each shower (9) are given in Table [2j 

To perform a background subtraction of the surviving cosmic ray events, an esti- 
mation of these backg round counts is made using the reflected-region background model 
(jAharonian et al.ll200ll ). Events within an angular distance 9 of the anticipated source lo- 
cation are considered ON events. Background measurements (OFF events) are taken from 
regions of the same size and at the same angular distance from the center of the FoV. For this 
analysis, a minimum of eight background regions was used. The excess number of events from 
the anticipated source location is found by subtracting the number of OFF events (scaled by 
the relative exposure, a) from th e ON events. S tatistical significances are calculated using 
a modified version of Eqn. 17 of iLi fc Ma (1983) to allow for varying number of off-source 
regions due to nearby star 60 Cygni (jAharonian et al.ll2004j ). More det ails about VERITAS , 
the calibration procedure and the analysis techniques can be found in lAcciari et al.l ( 120081 ). 



3. Results 

Analysis of the VERITAS data did not show a significant detection at the location of 
V407 Cyg. The results from both event reconstruction methods were used to c alcula te upper 



limits on the flux from V407 Cyg with the method described by iRolke et al.l (120051 ) and the 
assumption of a Gaussian-distributed background. 

The upper limits for V407 Cyg are calculated at the decorrelation energies of 1.8 TeV 
for the standard method and 1.6 TeV for the displacement method and assume that any 
emission takes the form of a power law with a photon index of -2.5. The decorrelation energy 
is the energy at which the dependence of the upper limit calculation on the assumed photon 
index is minimized. This energy is found by performing multiple upper limit calculations, 
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with different spectral indices, and determining the region where the resulting upper limit 
functions intersect. The energy threshold for the observations of V407 Cygni with VERITAS, 
defined as the maximum of the product of the assumed spectral shape and the effective area, 
is 1.2 TeV for both methods. 

The analysis results for V407 Cyg are presented in Table [3] and Figure [TJ In addition, 
results from observations of the Crab Nebula at similar zenith angles are presented in Table 
[3] and Figure HJ The efficiency of the displacement method for event reconstruction can be 
observed in the increase in both gamma-ray rate and significance for the Crab measurements. 
The increased sensitivity also results in a reduction of the upper limit for V407 Cyg compared 
to the standard method. 



Discussion 



The GeV detection of V407 Cyg provides evidence for previously unobserved gamma- 
ray emission from novae in white dwarf /red giant systems. Expanding shock waves have 
been known to accelerate particles to high energies, and gamma rays are observed from 
supernova remnants. The discovery by the Fermi-LAT team, however, suggests that the 
same phenomenon occurs in some novae, adding a new class of gamma-ray emitting objects. 
The lack of a significant detection in the VHE band suggests that either particles were not 
accelerated to sufficient energies to produce VHE photons during the V407 Cyg outburst or 
that VHE photons were produced, but then absorbed. 

The key to creating gamma rays is the acceleration of sufficiently energetic charged par- 
ticles. In the case of V407 Cyg, the expanding matter from the nova collides with the stellar 
wind from the red giant and causes a shock, which accelerates the particles near the shock to 
relativistic energies. A rough estimation for the maximum energy attainable by fi rst order 
Fermi acceleration of a particle at a shock can be found, following the discussion of iLongair 
( 1201 ll ). If B is the magnetic flux density where a shock proceeds and the shock travels with 
velocity U, the maximum energy of a particle with charge Ze is E max = J ZeBU 2 dt, where 
t is the time allowed for particle acceleration. This means that the highest attainable energy 
is proportional to the magnetic field in which the nova travels, the square of shock speed 
and the time for acceleration. The mean magnetic field in the shock can be estimated as 
B = [327T p(R)kT] 2 , where p(R) is the density of gas molecules wit h respect to the di stance 
from the center of the red giant and T is the wind temperature (lAbdo et al.l I201CH ) . We 
assume a wind te mperature of T=700K, corre s pondi n g to the temperature of t he dust enve- 
lope measured by lMunari. Margoni. fc Stagnil (119901 ) . lEspey fc Crowleyl (120081 ) were able to 
directly measure the temperature of the red giant wind in a similar symbiotic system, EG 
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And, and found that it can reach ~8000 K near the red giant. Using this wind temperature 
would increase the estimate for maximum particle energy by a factor of three. 



Orlando fc Drake! (1201 if ) carried out detailed hydrodynamic simulations of the V407 



Cyg nova with various gas distribution models and could accurately reproduce the X-ray 
light curve of V407 Cygni. Their model for the distribution of gas that best reproduced 
the light curve included what they call a "circumbinary density enhancement," a region of 
density exceeding the typical R~ 2 profile of the stellar wind in the binary system, and had a 
binary separation of 15.5 AU ( Orlando fc Drakelboil ). For th e temporal profile of the nova 
shock velocity, we used the equation that iMunari et al.l (120101 ) found from fits to the broad 
components of the Ha spectra they measured beginning at day 2.3 after the outburst (13 
March 2010) and thereafter: U = 4320-5440 logt+2635(logt) 2 -460(logt) 3 . For the velocity 
between day and day 2 of the outburst, we can assume two cases that bound the possible 
velocity profiles: i) the nova shell experienced free expansion at a constant velocity before 
day 2, with the assumption that the mass collected by the nova shell during this period was 
small (free expansion model) or ii) extrapolate the above equation for the velocity to times 
before day 2 (extrapolation model). We then find that at the start of VERITAS observations 
of V407 Cyg (day 9 of the outburst), E max ~ 1.4 TeV for the free expansion model and ~ 3.0 
TeV for the extrapolation model. This suggests that particles could have been accelerated 
to TeV energies by the time of the VERITAS observations. 

To check the importance of absorption, we calculated the opacity along the photon 
path for gamma rays generated at the shock front. Electron-positron pair p roduction via 
photon-photon collision is the dominant interaction ( iGould fc Schrederl 119671 ). We modeled 
the red giant spectrum as a blackbody with a temperature of 2500 K and a radiu s of 50 R Q , 
found the photon density as a function of position and ener gy following iNunea ( 2011 ). and 
used the cross section for the photon-photon collision from IGould fc Schrederl (119671 ) . We 
found that the opacity for TeV photons only becomes significant when the TeV emission 
region is located directly behind the red giant with the system viewed edge on. Though 
this case cann ot be ruled out, it is statisti c ally u nlikely. In addition, if the suggested orbital 
parameters of lMunari. Margoni. fc Stagnil (Il990l ) are accurate, the system is unlikely to have 
been in such an orientation at the time of the nova. 

The upper limits placed by VERITAS can put some restrictions on the gamma-ray 
emission mechanism in V407 Cyg . Two physical m odels of gamma-ray production at the 
shock- front have been suggested (lAbdo et al.ll2010l ). In the hadronic model, gamma rays 
are produced in the decay of 7r° particles generated by collisions of high energy protons 
accelerated in the shock. In the leptonic model, gamma rays are produced via inverse- 
Compton scattering of infrared photons emitted from the red giant on high energy electrons 
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accelerated in the shock. 



The electron threshold energy for production of a gamma-ray photon via inverse-Compton 



1 + 



scattering off the red giant photons can be estimated as: E > (E 7 /2) ^1 -| 
where E is the electron threshold energy, E 1 is the gamma-ray energy and e is the energy 
of the red giant photons. The electron threshold energy for a 1 TeV gamma ray scattering 
off 0.6 eV photons at the peak of the red giant spectrum is 1.1 TeV. Though the above 
calculation indicates that particles could reach TeV energies if continuously accelerated for 
the full nine days from the initial outburst to the start of the VERITAS observations, the 
inverse-Compton c ooling; time would be significantly less than a day (the time estimated by 
Abdo et al.l (12010[ ) for 5 GeV electrons), meaning electrons that are accelerated in the first 
few days of the outburst would not likely retain sufficient energy to produce VHE photons 
by the time of the VERITAS observations. VHE emission near the time of the VERITAS ob- 
servations would therefore require freshly accelerated particles, however, recently accelerated 
particles would likely not have enough time to reach TeV energies. If electrons did reach TeV 
energies, they would be approaching the Klein-Nishina regime, where the cross section for 
inverse-Compton scattering would be reduced, resulting in a longer cooling time. However, 
electrons with TeV energies would be well above the exponentia l cutof f, 2>.2^\ GeV, of the 
electron spectrum in the best-fit leptonic model of lAbdo et al.l (120101 ) . These two factors 
imply that the VERITAS upper limits place no new constraints on leptonic models. 

For hadronic models, the Fermi-LAT data alone provide relatively poor constraints on 
the extension of the proton spectrum to high energies. In particular, the cutoff energy (E cp ) 
is not well bounded from above if the spectral index is steep. The VERITAS data can be 
used to improve the con straints on t h e had ronic model parameters. To do so, we re-fit the 
hadronic model used by I Abdo et al.l (120 101 ) to the Fermi-LAT poin ts with the additi on of 
the VERITAS upper limit. Figure [3] shows the Fermi-LAT data (lAbdo et al.ll2010l ) and 
the VERITAS flux upper limit compared to the best fitting hadronic model. The gamma - 
ray spectrum is calculated via the method of iKamae. Karlsson. Mizuno. Abe, fc Koil ( 120061 ). 



-W p /E c 



assuming a cosmic proton spectrum of the following form: N p = N p fi(W p + m p c 2 )~ 
(protons GeV -1 ), where (W p , E cp ,m p are kinetic energy, cut-off energy, and mass of the 
proton and s p is the power law index). 

Figure H] shows a confidence region map for the parameters of the hadronic model using 
both the Fermi-LAT and VERITAS data. The gamma-ray spectrum was modeled as de- 
scribed above and fitted to the Fermi-LAT data by varying s p and E cp . The spectrum was 
then compared to the VERITAS upper limit, and a x 2 value for the VERITAS data point 
alone was calculated. Specifically, we calculated the model flux in the VERITAS energy 
band and compared this to the flux upper limit determined via the displacement method. 
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This x 2 value was then added to the \ 2 calculated for the Fermi-LAT data. The confidence 
levels were then calculated for the two parameters of interest, s p and E cp . As can be seen 
from Figure HI the VERITAS observations place greater restrictions on the model proton 
spectral index for high cutoff energies. The 90% confidence limits are E cp < 5 TeV (compa- 
rable to E max = 3 TeV, calculated above for the extrapolation model) and E cp > 0.01 TeV 
(much lower than E max ). It is possible that the peak energy of the particles pro duced by the 



shock could be reduced if the magnetic field is weaker than estimated above. iNelson et al. 



(120121 ) argue for a larger binary separation of 20-25 AU, based on the presence of lithium 
burning in the Mira, and evidence that the white dwarf in the system is massive. Using this 
larger separation distance would lead to weaker magnetic fields. The limits placed by the 
VERITAS observations are near the threshold for the observations, so it is also possible that 
simply not enough particles were accelerated to high enough energies to produce a significant 
detection by VERITAS. 
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Table 1:: Summary of the VERITAS data sets presented in this work. 



Target 


Period 


Useful Duration [min] 


Zenith Angle [°] 


V407 Cyg 


Mar. 19-26, 2010 


304 


50-66 


Crab Nebula 


Mar. 12-16, 2010 


203 


55-65 



This preprint was prepared with the A AS IATgX macros v5.2. 
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Table 2:: Selection criteria used for the VERITAS analysis. For an explanation of these 
parameters, see Section [2731 



Parameter 


Selection Criteria 


Image size 


> 400 digital counts (~ 75 photoelectrons) 


Mean Scaled Width 


0.05 < MSW < 1.15 


Mean Scaled Length 


0.05 < MSL < 1.3 


Height of Shower Maximum 


> 7 km 


e 


< 0.1° 



Table 3:: A summary of the VERITAS analysis results. Flux upper limits for V407 Cyg 
from VERITAS observations are calculated assuming a photon power-law index of -2.5 and 
taken at the decorrelation energy, minimizing the effect of the assumed spectral index. The 
displacement method produces both a higher significance for the Crab Nebula data and a 
more sensitive upper limit for V407 Cyg. 



Reconstruction method 


Standard 


Displacement 


Source 


Crab 


V407 Cyg 


Crab 


V407 Cyg 


Exposure [min] 


203 


304 


203 


304 


ON (Source) counts 


255 


91 


300 


76 


OFF (Background) counts 


744 


841 


351 


630 


a (See SecJ2.3p 


0.111 


0.125 


0.111 


0.125 


Significance 


15.7a 


0.5cr 


25.1a 


l.Ofj 


Rate [photons min -1 ] 


0.97 ±0.09 


0.02 ±0.04 


1.40 ±0.09 


0.03 ±0.03 


Energy threshold [TeV] 


1.5 


1.2 


1.7 


1.2 


Decorrelation energy [TeV] 




1.8 




1.6 


Flux upper limit at decorrelation 
energy (95% c.l.) [E 2 *dF /dE; erg 
cm -2 s _1 ] 




2.7 x 10~ 12 




2.3 x 10~ 12 
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V407 Cygni (standard) 



V407 Cygni (displacement) 
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(a) standard method 
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(b) displacement method 



Fig. 1. — : VERITAS significance maps of V407 Cyg reconstructed by the two methods. The 
white circle at the center indicates the location of the Fermi-LAT detection and the size of 
the source region used for the VERITAS analysis. 
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Fig. 2. — : VERITAS significance maps of the Crab Nebula reconstructed by the two meth- 
ods. 
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Fig. 3.—: Spec trum of V407 Cyg (FGL J2102+4542) measured by the Fermi-LAT 
( lAbdo et al.l 120101 ) and VERITAS upper limits. Vertical bars indicate la statistical er- 



ror, and arrows indicate 2a upper limit. The rightmost arrows show the 99% confidence 
level (3a ) VERITAS upper limit calculated using the displacement method (at 1.6 TeV) and 
the standard method (at 1.8 TeV) fo r event reconstruction (See section E)). The fitt ing curve 
was constructed with the method of iKamae. Karlsson. Mizuno. Abe, fc Koil (120061 ) with the 
parameters (s p , log(£ , cp )) = (2.15, 1.5) (See section 0J. 



(a) Fermi-LAT 



(b) Fermz-LAT+VERITAS 



Fig. 4. — : Confidence region maps for the hadronic model of gamma-ray production in V407 
Cyg constrained by Fermi-LAT data alone (a) and constrained by the Fermi-LAT data plus 
the VERITAS upper limit (b). The x-axis gives the spectral slope (s p ) and the y-axis is 
the logarithm of the cutoff energy (E cp in GeV). The numerical values correspond to the 
confidence level of the fit with the given values of the parameters. The 99% contour from 
Fermi-LAT data alone is shown as a dashed line in (b) for comparison. 



